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Crystal truncation rods calculated in the kinematical approximation are shown to quantitatively 
agree with the sum of the diffracted waves obtained in the two-beam dynamical calculations for 
different reflections along the rod. The choice and the number of these reflections are specified. The 
agreement extends down to at least ~ fO^^ of the peak intensity. For lower intensities, the accuracy 
of dynamical calculations is limited by truncation of the electron density at a mathematically planar 
surface, arising from the Fourier series expansion of the crystal polarizability. 
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The planar surface of a semi-infinite crystal gives rise 
to a diffraction pattern consisting of lines normal to the 
surface and passing through the bulk diffraction peaks. 
The intensity distribution along these lines, called crys- 
tal truncation rods (CTRs), contains information on the 
surface structurejii^ CTR calculations, commonly per- 
formed in the kinematical (single scattering) approxima- 
tion, are the basis of x-ray surface crystallography! '^I'^i^ '^ 
The high-intensity CTR regions close to the Bragg reflec- 
tions that require dynamical (multiple scattering) calcu- 
lations are not sensitive to surface structure and are usu- 
ally excluded from considerations in the surface structure 
determination studies. 

On the other hand, many studies of crystalline films 
and multilayers are based on dynamical diffraction theory 
and are mostly restricted to the close vicinity of the Bragg 
peaks J'-^'?' The substrate peak position is commonly used 
as a reference to study the film relaxation. The intensity 
of the substrate peak can also be used as a reference 
to obtain the film structure factor on an absolute scale 
to determine, e.g., ordering in sublattices.^'^ The thinner 
the film, the larger the part of the CTR that needs to 
be analyzed. Studies of very thin (a few atomic layers) 
films require the analysis of the whole CTR. The prob- 
lem of an accurate calculation of the CTR intensity in the 
whole wave vector range has been the subject of a num- 
ber of investigationsfiiii^ii^ii^ii^ii^iiiii^ii^ but a conclu- 
sive recipe has not been given yet. It is well established 
that the two-beam dynamical diffraction theory is very 
accurate in the vicinity of the respective Bragg refiec- 
tion but fails in the middle between two refiections. The 
resolution of the latter problem was attempted by appro- 
priate multibeam calculations i^^i^^d"^ Recently, Pavlov et 
al^ proposed to calculate the diffracted intensity in a 
wide angular range by summing up the amplitudes of 
the diffracted waves for different reflections, each ampli- 
tude being the solution of the corresponding two-beam 
diffraction problem. 

In the present paper, I show that the sum of the solu- 
tions of the two-beam diffraction problems for the reflec- 
tions along a CTR coincides with the kinematical for- 
mula. The number and the choice of reflections that 
have to be included in the summation are specified. Dy- 
namical and kinematical calculations disagree in the close 



vicinity of the Bragg reflections where the dynamical so- 
lution is correct while the kinematical one diverges, and 
in the regions where the diffracted intensity is less than 
~ 10^^ of the peak intensity. In this latter region, the 
dynamical calculations are less accurate. The error in 
the dynamical calculations originates from the Fourier 
expansion of the electron density over the reciprocal lat- 
tice vectors, so that it cannot be remedied by a more 
accurate multibeam dynamical theory. 
The analysis is based on the equality 
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It is derived below and the limits of its validity are es- 
tablished. Here hkl are the integer Miller indices of the 
refiections {h and k correspond to the directions in the 
surface plane and I is along the surface normal), L is 
the continuous coordinate along a crystal truncation rod 
(L = qa /2ti, where q is the z-component of the momen- 
tum transfer and a is the lattice spacing in z-direction), 
and FhkL — j f j ^'^vi'^T^iihx j + kyj -\- Lzj)/a\ is the 
structure factor of the unit cell (the sum is taken over all 
atoms in the unit cell, fj is the atomic scattering factor 
of the jth atom and Xj, yj, Zj are its coordinates). It will 
be shown below that Eq. ([T]) is valid if and only if the 
electron density obtained by the back Fourier transform 
of Fhki is entirely contained inside the unit cell. This 
latter requirement results in the cutting of the atomic 
electron densities at the surface. 

The left-hand side of Eq. ([1]) is proportional to the scat- 
tering amplitude of a semi-infinite crystal in the kinemat- 
ical approximation j^i^i'^i^ i^ The kinematical theory for 
semi-infinite crystals^ provides the remaining prefactors. 
The amplitude of the scattered wave is 
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Here A is the wavelength, is the classical radius of the 
electron, and 7out — sin<I>out, where ^out is the angle 
between the diffracted wave and the crystal surface. A 
similar factor for the incidence wave, 7in — sin<I>in, where 
^in is the angle between the incident beam and the sur- 
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face, will be used later. The amplitude of the incident 
wave is taken equal to 1. 

Each term of the sum on the right side of Eq. ([T]), after 
having been multiplied by the same factor Are/(a^7out), 
is the asymptotic expression of the two-beam dynam- 
ical diffraction theory for angular deviations from the 
Bragg peak that are much larger than the dynamical peak 
width. Within dynamical theory, the amplitude of the 
diffracted wave at large deviations from the Bragg peak 
for a given reflection H is given by Eu = XH/[27out(7in + 
Tout + V'h)]; where xh is the Fourier component of the 
polarizability and ipn = Hz/k (here k = 27r/A is the 
wave vector in vacuum). The position along the CTR 
L = (a/A)(7in + 7out) is expressed through the incidence 
and the exit angles, while the position of the Bragg peak 
on the CTR is I = —{a/X)ijjH- Substitution of these ex- 
pressions gives Eh = Arei^H/[27ra^7out(-^ — I)], which is 
equal to the corresponding term in Eq. ([1]) multiplied by 
the factor given above. Then, the sum of the amplitudes 
of diffracted waves calculated in the two-beam dynamical 
theory 

E''y'^{hkL)=Y,^Etlf (3) 

for all reflections hkl along a CTR hk coincides with the 
kinematical solution. For non-grazing incidence and exit, 
the resulting amplitude depends on the z-component of 
momentum transfer q = 2'KL/a, rather than the incident 
and the exit angles $in and $out separately. Pavlov et 
alr^ have arrived to an equivalent formula by solving the 
Takagi-Taupin equations. 

In the close vicinity of each reflection hkl^ for angular 
deviations comparable to the Darwin width, the ampli- 
tudes of the other reflections are much smaller (on the 
order of the polarizability x ^ 10~^) than that of the 
actual reflection, and Eq. ([3]) can be applied, keeping in 
mind that all other terms can be neglected compared to 
^hkt- larger angular deviations, Eq. ([3]) is the flrst 
order (in x) approximate solution of the multibeam dy- 
namical diffraction problem that includes all Bragg re- 
flections along a given CTR. Such an approximation is 
valid since at most one Bragg reflection occurs at once 
at the diffraction condition when moving along the CTR. 
Hence, Eq. ^ describes the intensity distribution along 
a CTR around Bragg peaks and far from them, down to 
intensities ~ x of the peak intensity, which covers almost 
the whole CTR. Eq. ([3]) fails at the order of x^- Before 
ascertaining the problems arising at such low intensities, 
let us consider the CTR calculations based on Eq. ([3]). 

I have implemented the dynamical diffraction equa- 
tions in the formulation of Stepanov et al.^^ Since Eq. 
([5]) involves the summation of the scattering amplitudes, 
a proper treatment of the phases of the Fourier compo- 
nents Xhki of the polarizability is essential. The structure 
factor calculations require the origin of the unit cell to 
be chosen once for all reflections. This is not so in the 
programs by Sergey Stepanov^^i that choose the origin 
separately for each reflection to minimize the phase of 



Xhki- Such a choice is correct as long as a single re- 
flection is involved in the calculations. When working 
with the amplitudes of several reflections, Eq. this 
is not appropriate. Also, for layered structures, the in- 
terference pattern depends on the relative displacements 
of the crystal lattices of different layers, which requires 
a common origin for the calculation of all Xhki'sr^ The 
requirement of a common origin is usual in surface x-ray 
structure analysis but not appreciated in the dynamical 
diffraction calculations. 

Equally important, but probably less evident, is the 
choice of the atomic positions in the unit cell with re- 
spect to the surface. The surface is commonly taken at 
z = 0, and one of the atoms is placed at the origin. If 
so, the surface cuts half of this atom (and other atoms at 
the same level), which leads to wrong CTR intensities. 
For example, with the surface at z = 0, the Ga and As 
atoms in the bulk GaAs unit cell cannot be placed at the 
depths 0, 1/4, 1/2, 3/4. Rather, they have to be put at 
the depths 1/8, 3/8, 5/8, 7/8. Only this choice leads to 
an agreement between dynamical and kinematical calcu- 
lations. Note that this choice of the atomic positions is 
irrelevant to the kinematical formula ([U, since the kine- 
matical calculation, based on the summation over atoms 
and unit cells, contains only one structure factor F^kL 
(with the continuous parameter L), the phase of which 
is lost when calculating the intensities. In contrast, the 
relative phases of the structure factors Fhki (with integer 
Vs) depend on the choice of atomic positions. 

The choice of the reflections that have to be included 
in the sum ^ follows from Eq. ([T]). Since L starts from 
zero, positive and negative Vs are equally important. 
Hence, the sum should contain the solutions of the two- 
beam dynamical problems in both reflection (the Bragg 
case) and transmission (the Laue case). Although the 
dynamical diffraction equations are formulated in Ref . 
for the reflection case, they do not require any change to 
include the transmission case, since the diffracted wave 
is anyway directed back into vacuum. In other words, 
the reflections that formally correspond to the Laue case 
(the diffracted wave propagates into the crystal when the 
Bragg condition is satisfied) become the Bragg case re- 
flections because the deviation from the Bragg condition 
is so large that the diffracted wave propagates into the 
vacuum. The number of equations to be included in 
the sum ([1]) or ([3]) depends on the decay of the struc- 
ture factors Ffiki with increasing I. The atomic scatter- 
ing factors /(s) are decreasing functions of s = sin6'/A, 
where 9 is the Bragg angle of the respective reflection, 
with s ~ 2 as a characteristic scale for this decay. 
Then, for reflections with I ^ h,k the lattice spacing is 
d — a/l and, from Bragg's law s = l/2(i, reflections up to 
^max ~ 4a have to be included (here a is in angstroms). 
Hence, a reasonable estimate is Zmax ~ 25. In fact, the 
calculations below show that l-ma^ ~ 10 is sufficient for 
practical calculations. 

Figure [Ua) compares kinematical and dynamical cal- 
culations of the CTR IIL from a hypothetical bulk- 
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FIG. 1: Crystal truncation rods IIL from a bulk-terminated GaAs(OOl) crystal at (a) non-grazing incidence angle "l>o = 5° 
and (b) grazing incidence angles <l>o = 0.3° and $o ~ 0.2°, above and below the critical angle Qc — 0.24°. (c) CTR IIL from 
a 56.53 A thick FeaSi film on GaAs(OOl) at an incidence angle "l>o = 0.5°. The two right panels magnify the curves of the left 
panel at the intensity minimum (bulk- forbidden reflection 112) and close to the maximum (bulk refiection 111). Dynamical, 
kinematical, and distorted-wave Born approximations are compared. 



terminated GaAs(OOl) crystal. The reflectivities R — 
(7out/7in) \E{hkL)\ are plotted. A non-grazing inci- 
dence angle $o — 5° is chosen to ensure kinematical 
diffraction conditions. In the left panel, presenting the 
whole CTR, the two calculated curves are almost indis- 
tinguishable. At the 111 bulk reflection (the right panel) 
the kinematical curve is slightly shifted with respect to 
the dynamical one because of refraction at the surface, 
and diverges at the Bragg position. A small but impor- 
tant discrepancy between the curves is revealed in the 
middle panel that enlarges the region of the intensity 
minimum, close to the bulk-forbidden reflection 112. Dy- 
namical calculations are performed with different num- 
bers of terms in the sum ([3]). One can see from the plot 
that 10 terms (i.e., the sum of two-beam dynamical cal- 
culations for the reflections 111 with odd I in the range 
— 9 < / < 9) are sufficient for the convergence of the se- 
ries. However, the sum ^ does not converge to the kine- 
matical solution even if the number of the involved reflec- 



tions is increased further. The origin of the disagreement 
is in the Fourier expansion of the electron density over 
the reciprocal lattice vectors hkl. This expansion is cer- 
tainly correct for an inflnite bulk crystal that consists of 
a periodic repetition of the unit cells. However, if a unit 
cell is cut out of the crystal, its electron density does not 
coincide with the sum of electron densities of the atoms 
that belong to this unit cell. Figure [2] illustrates this 
statement. Here the atoms of the GaAs bulk structure 
are represented by balls with radii equal to their covalent 
radii. The broken line selects one unit cell. Its position is 
chosen to minimize the overlap of the electron densities 
belonging to the atoms of different unit cells. Still, parts 
of the electron density of atoms of the chosen unit cell are 
out of that unit cell, and parts of atoms from neighboring 
unit cells fall into the chosen unit cell. These do not af- 
fect the bulk structure factor calculations, since the unit 
cells are periodically repeated. However, the dynamical 
Bragg diffraction from a scmi-inflnite crystal implies a 
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FIG. 2: View of a truncated GaAs bulk structure in the [100] 
direction. Atoms are represented by balls with their radii 
equal to the covalent radius of the respective element. A unit 
cell marked by broken lines is chosen to minimize the overlap 
of the electron densities from atoms of different unit cells. 
The dotted straight line represents a cut by a mathematically 
planar (001) surface. The solid line is a physical (001) surface 
terminated by atoms. 



truncation of the sum over the unit cells, which is equiv- 
alent to a rigid truncation of the electron density by a 
plane (dotted line in Fig. ^ . Such a cut of the electron 
density of an infinite crystal by a plane removes parts 
of the top atoms and artificially adds the parts of the 
atoms belonging to the next layer. In contrast, the phys- 
ical surface is rough on that scale since it contains the 
entire densities of the atoms of the top layer, as shown 
by the thick line in Fig. [21 

The difference in these electron densities does not play 
a role in the vicinity of Bragg peaks, where many atoms 
contribute to diffraction, but becomes essential between 
Bragg peaks, where the surface atoms provide the main 
contribution to the diffracted wave. To provide an ad- 
ditional proof that this effect is the only source of dis- 
crepancy between the dynamical and kinematical calcu- 
lations, I have repeated the same calculations for artificial 
atoms two times smaller in size placed in the same lat- 
tice, Fig. [3l Such a crystal is not, in fact, physical, since 
the long-range attractive forces between atoms are not 
compensated by the short-range repulsion when the dis- 
tances between atoms exceed the sum of their covalent 
radii by a factor of 2. The atomic scattering factors for 
the artificial atoms Xa and Xb are chosen to be Gaussians 
with widths approximately two times larger than the /(s) 
functions for Ga and As atoms. Fig. El^b). The conver- 
gence of sum ([T]) or ^ is fairly slow, Fig.[3^c). The sum 
of dynamical amplitudes converge to the distorted-wave 
Born approximation, described below, when 20 terms are 
included. Each term is a solution of the two-beam dy- 
namical diffraction problem for a reflection IIZ with odd 
I from -17 to +21. 

This result discourages other improvements of the dy- 
namical theory, in particular a complete solution of the 
multibeam diffraction problem, that could improve the 
accuracy compared to Eq. ([3]). As long as a Fourier ex- 



pansion over reciprocal lattice vectors is employed, the 
cut of the electron density shown in Fig. [5] persists and 
does not provide an accuracy that is better than given 
by Eq. ([3]), i.e., an accuracy of the order of ~ 10~^ with 
respect to the peak intensity. Fortunately, lower intensi- 
ties quite rarely arise in experimental CTR studies and 
are additionally masked by surface reconstruction and 
roughness. 

Equation ([3]) can be directly applied to grazing inci- 
dence diffraction with the same reasoning as above: in 
the vicinity of a Bragg peak, corrections due to other re- 
flections are negligible, and in the remaining part of a 
CTR the sum ([3]) is the first order (over x) perturba- 
tion solution of the multibeam diffraction problem. The 
kinematical approximation is extended, for the grazing 
incidence diffraction the distorted- wave Born ap- 

proximation (DWBA).^i22,23^ rpj^g DWBA formulation 
in terms of the reciprocity theorem in electrodynamics is 
most straightforward3»2i In the zcroth order, the scat- 
tering problem is solved for a uniform medium having the 
same polarizability xo a-s the crystal under investigation. 
The scattering problem is solved twice, with the waves in- 
cident on the surface under the incidence angles and 
'i'out, respectively. The solution of each problem in the 
medium consists of two plane waves with the amplitudes 
Di [i = 1,2), corresponding to the transmitted and the 
specularly reflected waves in vacuum. A convenient way 
to find these solutions is to reduce the 4x4 matrices of 
the dynamical diffraction problemSfl to 2x2 matrices. In 
this way, layered structures with different xo's can easily 
be treated. Then, the kinematical solution is replaced by 

i?°'^^^($i„,<fout) = J2 DTD^'^'E'^^ihkL,,). (4) 

The superscripts "in" and "out" distinguish the respec- 
tive zeroth order solutions. The parameters = 
(a/A)(u-" 4- u™') are obtained from the parameters Ui — 
kzi/K describing the wave vectors in the medium (here 
ki are the wave vectors of the waves inside the medium). 
Complex parameters Lij substitute the real parameter L 
in Eq. ^. 

Figure [IJb) compares dynamical and DWBA calcula- 
tions of the same CTR IIL for the incidence angles below 
and above the critical angle ac = 0.24°. A good overall 
agreement is evident from the left panel. The discrepancy 
at the intensity minimum (middle panel) has the same 
nature as above. The right panels of Fig. [2a,b) show 
that the DWBA and the dynamical calculation agree ev- 
erywhere except within the Darwin width of the Bragg 
peak. If the incidence angle is smaller than the critical 
angle, the results of the two calculations coincide in that 
region as well. 

The extension to layered structures is straightforward. 
The diffracted beam amplitudes obtained in the two- 
beam dynamical calculations^^ are summed up accord- 
ing to Eq. ([3]). The kinematical amplitude ^ is sup- 
plemented by an additional term arising from the fi- 
nite sum over the layer's unit cells, f ^^.^ [exp(27rzA^L) — 
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FIG. 3: (a) An artificial crystal XaXb with tiie same structure and the lattice parameter as GaAs but with atomic sizes 
that are two times smaller. The unit cell (broken lines) and a planar surface (dotted line) can be cut between atoms, (b) 
Atomic scattering factors attributed to atoms Xa and Xb. The functions f{smd/X) are chosen to be Gaussians with widths 
approximately two times larger than those of Ga and As. (c) The CTR IIL from the XaXb semi- infinite crystal in the vicinity 
of L « 2 calculated within distorted- wave Born approximation (DWBA) and dynamically, by Eq. (|3]), with different numbers 
of involved refiections. As the number of reflections is increased, the dynamical curves converge to the DWBA results. The 
incidence angle is <l>o = 0.5°. 



l]/[exp(27rjL) — 1], where -F^j,^ is the structure factor of 
the layer unit cell, and Na is the layer thickness. Figure 
[Ijc) compares dynamical and DWBA calculations for a 
FesSi film on GaAs(OOl). FeaSi possesses a cubic unit cell 
with the same lattice spacing as GaAs, so that there is no 
mismatch between the two lattices.^*' The incidence angle 
is taken $0 = 0.5°, so that the DWBA corrections to the 
kinematical formulas are essential. The film is = 10 
unit cells thick {Na = 56.53 A). A perfect agreement 
between the two curves is evident. Even at the inten- 
sity minimum (middle panel), the discrepancy between 
the two calculations is absent. At the bulk reflection 111 
(right panel), the DWBA solution diverges and cannot 
be used within the Darwin width, while the dynamical 
solution is correct. The dynamical calculations in Figs. 
[IJb,c) were performed by the summation of solutions of 
32 two-beam diffraction problems (11/ refiections with 
odd I in the range —31 < / < 31 were calculated). 

The CTR IIL analyzed in Fig. [1] shows a larger dis- 
crepancy between dynamical and kinematical calcula- 
tions than other rods from GaAs(OOl) surface. Figures 
|4lja,b) present two other rods. The difference between 
the two calculations diminishes for the 13i rod and al- 
most disappears for the 22L rod. Further calculations 
show that the discrepancy decreases for AniBy crystals 
consisting of elements of two different periods in the pe- 
riodic table [e.g., GaSb, see Fig. (He)] and increases for 
lighter elements from the same period [e.g., AlP, see Fig. 
Sljd)]. The discrepancy remains quite small on the scale 
of the whole CTR and hardly can have practical implica- 
tions. Elemental semiconductors show more sharp min- 
ima at the bulk-forbidden reflection 112, see Figs. [ll^e,f). 
Similarly to the case of compound semiconductors, sili- 
con as a lighter element shows a larger, as compared to 
germanium, discrepancy between dynamical and DWBA 
calculations, see Fig.HJ^f). 

Equation ([T]) remains to be derived. Consider an arbi- 
trary complex one-dimensional function p{z) that is equal 



to zero outside a finite interval [0, a]. The Fourier integral 
of the function p{z) therefore involves integration over a 
finite interval: 
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On the other hand, the function p(z) can be periodically 
repeated along the z-axis. 
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and the periodic function Poc{z) can be expanded into 
the Fourier series: 



,{z)^ ^ F{qi)exp{-iqia), 



(7) 



where qi = 27:1 /a. Let Q.{z) be a function equal to one on 
the interval [0, a] and equal to zero outside this interval. 
Then, identically 



p{z) = Poo{z)^{z)- 



(8) 



The Fourier integral of the left-hand side of this equation 
is equal to F{q). The Fourier integral of the right-hand 
side is calculated by performing the integration for each 
term of the sum: 

F{<l) = / exp[i{q - qi)a]dz 
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It is taken into account that exp(iqia) — 1 since qi = 
27:1 /a with an integer I, so that qia is a multiple of 
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FIG. 4: Crystal truncation rods from GaAs (a,b), other AmBy compounds (c,d), and elemental semiconductor (e,f) bulk- 
truncated (001) surfaces. Dynamical (full lines) and DWBA (broken lines) calculations are compared. The incidence angle is 
$0 = 0.5°. 



27r. The function p{z) can now be identified with the 
Fourier component phk{z) in the expansion of the elec- 
tron density of a one unit cell thick crystalline slab over 
the wave vectors of its two-dimensional reciprocal lat- 
tice. Its Fourier transform F{q) is equal to FhkL with 
L = qajl-K, and F{(ii) are equal to Fh]^i. Equation Q 
reduces to Eq. ([T]). The crucial point in the derivation 
is contained in Eq. ([8]), which explicitly requires that the 
function p{z) is restricted to the interval [0, a]. 

Hence, the structure factors in Eq. ^ correspond to 
the electron density constrained to a unit cell cut out 
of the crystal, rather than to the electron density of the 
atoms whose centers are inside this unit cell. Parts of the 
electron density distributions of the atoms belonging to 
the chosen unit cell are cut away, while parts of the atoms 
from the surrounding unit cells occur in the chosen unit 
cell. Fig. O The problem can be avoided in rare cases 
of layered crystals, e.g. graphite, but is enhanced for 
higher- index surfaces. This effect limits the applicability 
of the Fourier series expansion of the electron density for 
the solution of the diffraction problems. The solution of 
a multibeam diffraction problem instead of Eq. ([3]), as 
well as corrections to the dynamical equations that are 
omitted in Ref . [20l . provide improvements of the order of 
X^, smaller than the effect of rigid truncation. That is 
why these corrections are not included in the calculations 
here. However, only the regions of extremely low inten- 
sity are affected. Intensities in these regions are sensitive 



to surface reconstruction and roughness. Hence, further 
improvements of the dynamical theory seem of limited 
practical impact nowadays. 

In conclusion, it is shown that the kinematical calcu- 
lation (or, at small incidence angles, the distorted-wave 
Born approximation) quantitatively agrees with the sum 
of the diffracted beam amplitudes obtained in the two- 
beam dynamical calculation. The number of Bragg re- 
flections that have to be included in the dynamical cal- 
culations is estimated from the angular dependence of 
the atomic scattering factors and amounts to some tens 
of reflections. Both transmission (Laue) and reflection 
(Bragg) cases have to be included. The reference unit 
cell should be chosen to minimize (since it cannot be 
completely excluded) the cutting of the electron densities 
of the top atoms by the surface. The agreement between 
kinematical and dynamical calculations can be lost in two 
regions. In the Darwin width regions near Bragg reflec- 
tions the dynamical theory provides correct intensities 
while the kinematical theory diverges. In the regions of 
very low intensity (below ^ 10~^ of the peak intensity), 
dynamical theory may fail because of the electron den- 
sity truncation by a mathematically flat plane instead of 
the physical surface while the kinematical theory remains 
applicable. 

This work has been inspired by discussions at the 8th 
Conference on High Resolution X-Ray Diffraction and 
Imaging (XTOP 2006). 
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